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ABSTRACT 



Aims. We present an analysis of physical conditions in planetary nebulae (PNe) in terms of collisionally-excited line 
(CEL) and optical-recombination line (ORL) profiles. We aim to investigate whether line profiles could be used to study 
the long-standing CEL/ORL abundance-discrepancy problem in nebular astrophysics. 

Methods. Using ID photoionization models and their assumed velocity fields, we simulate the line profiles of various 
ionic species. We attempt to use our model to account for the observed CEL and ORL profiles. As a case study we 
present a detailed study of line profiles of the low-excitation planetary nebula (PN) IC 418. 

Results. Our results show that the profiles of classical temperature and density diagnostic lines, such as [O ill] 
AA4363,5007, [S ll] AA6716,6731, and [Ar iv] AA4711,4740, provide a powerful tool to study nebular temperature and 
density variations. The method enables the CEL/ORL abundance-discrepancy problem to be studied more rigorously 
than before. A pure photoionization model of a chemically-homogeneous nebula seems to explain the observed disagree- 
ments in the profiles for the [O ill] A4363 and the A5007 lines, but cannot account for the differences between the [O ill] 
CELs and the O ll ORLs. We also investigate the temperature and density variations in the velocity space of a sample 
of PNe, which are found to be insignificant. 
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1. Introduction 

As descendants of AGB (asymptotic giant branch) stars, 
planetary nebulae (PNe) serve as crucial probes of the 
chemical evolution of galaxies and stellar nucleosynthetic 
processes. However, the physical conditions inside PNe 
are not yet completely understood. An important prob- 
lem in nebular astrophysics is that the heavy-element abun- 
dances measured from collisionally-excited lines (CELs) are 
systematically lower than those determined from optical- 
recombination lines (ORLs) (see Liu | 2006al for a recent 
review). The most extreme known case is that of Hf 2-2 
which has an abundance discrepancy- factor (ADF) of about 
70 (Liu et al. I2006p . Possible explanations include electron 
temperature, electron density, and chemical-abundance in- 
homogeneities (Peimbert 1967; Viegas & Clegg 1994; Liu 
et al. I2000p . implying that PNe have complex physical con- 
ditions. 

Many efforts have been made to develop new meth- 
ods to probe nebular physical conditions. The traditional 
methods of using strong forbidden lines as plasma diagnos- 
tics suffer some disadvantages. Because of their low critical 
densities (^ 10^ cm~'^), CELs are strongly suppressed by 
collisional de-excitation in high-density regions, and can- 
not be used to probe the physical conditions inside dense 
clumps within the nebulae. Given that the fluxes of CELs 
have an exponential dependence on electron temperature, 
they also cannot be used to probe extremely-cold nebu- 
lar regions where strong ORLs may arise. Improvements 
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in the quality of observations and increasingly-accurate 
atomic data are allowing us to use ORLs as plasma di- 
agnostics of PNe. Zhang et al. (|2004p presented a method 
to obtain nebular electron temperature and density simul- 
taneously using hydrogen-recombination spectra. Applying 
the method to 48 Galactic PNe, they found that the den- 
sities derived from hydrogen-recombination spectra were 
generally higher than those derived from forbidden line ra- 
tios and that the temperatures deduced from hydrogen- 
recombination spectra, To(H i), were systematically lower 
than those derived from the [O ill] nebular-auroral line- 
ratios, To([0 III]), suggesting that temperature and den- 
sity variations are generally present in PNe. Zhang et al. 
P005a( I2005bl) used He i recombination line-ratios to de- 
termine nebular electron temperatures. To (He i), and found 
that re(He i) < Te(H i), supporting the hypothesis that 
some cold hydrogen-deficient clumps are embedded in dif- 
fuse gas with a "normal" chemical composition, as proposed 
by Liu et al. ()2000p to explain the CEL/ORL abundance 
discrepancies. Some researchers presented nebular temper- 
atures and densities derived from faint O ii recombination- 
lines (Liu l2006bj Peimbert & Peimbert I2005|, and references 
therein) . Their results suggest that O ii recombination- lines 
may arise from colder and denser regions compared to the 
diffuse nebular gas. 

Line profiles of PNe are of primary importance to inves- 
tigate nebular physical conditions. To understand nebular 
dynamical evolution and geometry, Gesicki et al. (200]3]) 
and Morisset & Stasihska (|2006aP studied velocity profiles 
of a few strong emission lines, in terms of ID and 3D pho- 
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toionization models. A catalog of emission line profiles for 
PNe with various morphologies was presented by Morisset 
& Stasinska (2006b, 2008). Emission line profiles can also 
provide an opportunity to test whether CELs and ORLs 
originate in different zones. To complete such a test, nebular 
spectroscopic data of high resolution and signal-to- noise are 
required. Ruiz et al. ()2003l) and Peimbert et al. (2004) anal- 
ysed echelle spectra of the PNe NGC 5307 and NGC 5315. 
They found that O ii ORLs and [O ill] CELs had simi- 
lar widths, implying that both lines originated in the same 
volume and disfavouring the scenario of chemical inho- 
mogeneities. However, high spectral-resolution observations 
of two PNe NGC 6153 and NGC 7009 which have large 
CEL/ORL abundance discrepancies demonstrated that O ii 
ORLs are significantly narrower than [O ill] CELs (Barlow 
et al. [20061 Liu l2006ap . implying that the ORLs and CELs 
originate in different zones. The spectra of NGC 6153 ob- 
tained by Barlow et al. (t2006|} show that the [O ill] A4363 
and A5007 lines have remarkably-different profiles, indicat- 
ing the presence of large temperature variations within the 
nebula. 

In this paper, we simulate nebular emission line pro- 
files using photoionization models. Our main objective is 
to investigate whether line profiles could be used to study 
nebular physical conditions and the problem of ORL/CEL 
abundance-discrepancies in PNe. This paper is organized 
as follows. Section 2 presents our method. In Sect. 3, we 
use line profiles to study temperature and density variation 
within nebulae and compare with observations. Discussions 
and conclusions are presented in Sect. 4. 

2. Method 

The basic method is similar to that used by Gesicki et al. 
(|2003|) and Morisset & Stasinska (2006a). We first construct 
photoionization models using the CLOUDY code (Ferland 
et al. I1998p . A black-body spectral energy distribution is 
adopted for the ionizing stars. The nebulae are assumed 
to be a spherical shell with a given radial density distribu- 
tion. CLOUDY calculates the emissivity distribution of the 
CELs and ORLs for different ionic species and the ther- 
mal structure of nebulae. Using known slit aperture and 
positions and based on an assumed radial-velocity field, we 
deduce the integrating line profiles by summing along the 
line of sight. For our modellings, we assume that the slit 
width is far narrower than the nebular radius since high- 
resolution spectra are, in all cases, essential for the study. 
The predicted line profiles are then derived by convolving 
the broadening caused by temperature, seeing conditions, 
and instrument. 

Turbulence is not taken into account in our modelling. 
Gesicki & Zijlstra (2003) studied the expansion-velocity 
fields of a number of PNe and found that turbulence might 
be significant in Wolf-Rayet PNe. Based on a 3D modelling, 
Morisset & Stasinska (|2006a|) . however, argued that no tur- 
bulence is required to explain the observed line profiles if 
one considers a nebular-model departure from the assump- 
tion of spherical symmetry. It is therefore safe to assume 
that the effects of turbulence, on a normal PN, can be ne- 
glected. 

Figure [T] provides an illustration of our modelling of 
a medium-excitation PN. The velocity field is assumed to 
have a positive gradient across the nebula. An inspection 
of Fig. [T] shows that the high-ionization lines are narrower 
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Fig. 1. Modelled [O i] (first row), [O ll] (second row), and 
[O III] (third row) fine profiles for a nebula observed by a 
long-slit located across the nebular centre (first column), 
half of the nebular shell (second column), and the nebular 
edge (third column). The peak flux of each line is normal- 
ized to unity. Solid, dashed, and dotted lines correspond to 
three different velocity fields shown in the lower panel. 



than the low- ionization lines. This is because the high- 
ionized species are mainly located within the internal re- 
gions where the radial velocity is lower. The profiles of emis- 
sion lines arising from species with different ionization po- 
tentials have different dependencies on the velocity fields, 
and provide constraints of the nebular dynamics. Gesicki 
et al. (fT996l fT998|) and Gesicki & Zijlstra (20001 120031) in- 
vestigated velocity fields for a sample of PNe by fitting 
the strong H i, [N ii], and [O ill] line profiles. They found 
that the expansion velocities generally increase outwards, 
consistent with the predictions of hydrodynamical calcula- 
tions. Deep spectroscopic observations enable weaker lines 
including CELs and ORLs which cover a wider range of 
ionization potentials to be detected, and tighter constrains 
on the determination of velocity fields to be made. In the 
deep echelle spectra of the young PN IC 418, Sharpee et 
al. (|2004p detected profiles of a number of CELs and ORLs 
and found a negative correlation between line widths and 
ionization potentials. This is an implication that nebu- 
lar expansion- velocities increase outwards. Throughout this 
paper, a monotonously-increasing velocity field is adopted. 

For the model shown in Fig. [l] the [O ii] line profiles 
are more sensitive to the velocity field than the [O i] and 
[O III] line profiles because the 0+ species are the dominant 
ionization-state for the modelled medium-excitation neb- 
ula. Figure [H also shows that the line widths decrease when 
the slit is placed towards the nebular edge. This tendency 
is particularly evident for low-ionization lines originating in 
external nebular regions. For a given expanding shell, this 
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Fig. 2. Modelled results illustrating the effect of tempera- 
ture variations on line profiles: a) the predicted line profiles 
of [O III] A4363,[0 III] A5007, and O ii A4152; b) the as- 
sumed velocity field; c) the radial temperature distribution; 
d) the radial 0^+/0 abundance-ratio distribution. 



tendency is present because the velocity projected along 
the line of sight decreases from nebular centre to edge, and 
velocity broadening becomes less important with distance 
from the nebular centre. At the outer edge of nebulae, ther- 
mal broadening is dominant. Therefore, the emission lines 
originating in the outer edge of nebula have a Gaussian 
profile which is essentially independent of expansion veloc- 
ity and only depends on electron temperature and atomic 
mass, as shown in Fig. [TJ 

3. Analysis 

3.1. Emission line profiles as a probe of temperature 
variations 

3.1.1. Modellings 

Two scenarios were proposed to account for the ORL/CEL 
abundance discrepancies: 1) temperature fluctuations 
within the nebula fPeimbert 11967^ : 2) the two-component 
nebular model with extremely-cold hydrogen-deficient in- 
clusions embedded in the diffuse nebula (Liu et al. I2000p . 
The investigation of nebular temperature structure is vital 
for the solution of this problem. To quantitatively char- 
acterize the temperature variations within PNe, Peimbert 
(|1967p introduced the average temperature Tq and the 
mean square temperature fluctuation parameter , defined 

by 

jTMhdV 



and 

2 _ /(Te - To)^N,N,dV 

TSJNMV ' ^' 

respectively, where iVj is the number density of the ionic 
species used to determine the electron temperature. For a 
nebula with small temperature fluctuations, Tq and can 
be deduced by comparing Te([0 ill]) and Te(H i). By study- 
ing a large number of PNe, Zhang et al. (|2004p found, how- 
ever, that the average t^-value of the nebulae was higher 
than predicted by photoionization models. This is proba- 
bly evidence of the presence of extremely-cold 1000 K) 
clumps in which ORLs are significantly enhanced and CELs 
are barely excited. To study nebular temperature fluctua- 
tions, Rubin et al. (|2002p obtained a map of the Tc([0 ill]) 
of the PN NGC 7009 on the basis of high spatial-resolution 
Hubble Space Telescope images. The spatial variation of 
re([0 III]) in a number of PNe was also studied by Krabbe 
et al. (|2DD5)) . They found that the value of To([0 ill]) re- 
mains almost uniform across the PNe. Nevertheless, their 
results cannot completely rule out the scenario of temper- 
ature fluctuations as the cause of ORL/CEL discrepancies 
since their measurements of ^^.([O III]) are an average along 
each line of sight and provide only a lower limit to i^. A 
comparison of line profiles may provide an unique way to es- 
timate the temperature fluctuations along the line of sight. 

We study the effect of temperature variations across the 
nebula on the profiles of the [O ill] A4363 and A5007 lines. 
Figure [2] shows the modelled results of a high-excitation PN 
in which O^"*" is the dominant ionization state. This model 
assumes that the hydrogen density in the nebula is n(R) — 
10^ cm~'^ and that the central star has an effective tem- 
perature of 150,000K and a luminosity of 4 x lO'^^ergs"^. 
As shown in Fig. [2l we assume that the velocity field fol- 
lows a Hubble law. The slit is assumed to be midway be- 
tween the central star and outer edge of the nebula. For 
the models used throughout this paper, the average chemi- 
cal composition of the Galactic PNc (Kingsburgh & Barlow 
I1994P is assumed. As shown in Fig. [2l the predicted elec- 
tron temperature first decreases along the radius and then 
slightly increases at the outer edge of the nebula. The low- 
temperature regions are tightly associated with the 0^+/0 
abundance ratio. This is a natural consequence of the rela- 
tively strong [O III] nebular lines that provide the dominant 
means of cooling in this PN. 

Figure [2] shows that the variations of [O ill] A4363 
and A5007 line profiles are very sensitive to nebular ther- 
mal structure, and thus provide an opportunity to probe 
nebular temperature variations. In the low-velocity regions 
where the electron temperature is higher than the average 
value, the [O iii] A4363 line has a higher normalized inten- 
sity than the [O ill] A5007 line. This is because the auro- 
ral forbidden line has a higher excitation energy, providing 
higher weighting to high temperature regions. 

Based on high spectral resolution (i? = 150, 000) obser- 
vations of the PN NGC 6153, Barlow et al. ({2006)) found 
that the [O ill] A4363 and A5007 lines have different pro- 
files, that is, the gap between the two peaks of the A5007 
line is filled in the case of the A4363 line. This can be ex- 
plained in terms of our modelling, as shown in Fig. [2] It 
is interesting that the value 0.01) required to match 
the observed [O ill] AA4363, 5007 line profiles is much lower 
than that derived from the H i Balmer Jump and [O ill] 
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Table 1. The values" of Tq ^ and tl. 
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Fig. 3. Electron-temperature distribution in the velocity 
space of PNe, IC 2501, IC 4191, and NGC 2440. The dashed 
lines represent the values of Tq^^. 

lines (t^ = 0.045; Liu et al. I2000p . This appears to sup- 
port the idea that the two-abundance nebular-component 
model accurately describes the observations, to discount the 
theory that temperature fluctuations are the cause of the 
CEL/ORL abundance dichotomy. These conclusions are, 
however, arguable since the analysis of [O ill] AA4363, 5007 
line profiles can provide only a lower limit to the tempera- 
ture variations along the line of sight (see below). 

3.1.2. Temperature variations in the velocity space 

To study the temperature variations, we define the average 
temperature Tq^i, and the mean square temperature fluctu- 
ation parameter in velocity space by 



and 



/ T,{v)h{Y^(3)dv 



J[n{v)-ToJ^I,{mdv 



(3) 



(4) 



where v is the velocity along the line of sight, and Tc{v) 
and ly (up) are the electron temperature deduced from the 
forbidden line intensity ratio and the intensity of H/3 (cx 
NpNc) at a given velocity, respectively. The integrations 



Object 


IC 2501 


IC 4191 


NGC 2440 


To (K) 


7700 


8600 


12600 


0.053 


0.040 


0.077 


To.„ (K) 


11500 


10200 


14800 


0.002 


0.001 


0.011 



For the low-excitation PN, IC 2501, Tq^v and are deduced 
using the [N ll] AA5754, 6678 lines; for the high-excitation 
PNe, IC 4191 and NGC 2440, To,„ and tl are deduced using 
the [O III] AA4363, 5007 lines. 



in the two equations are over velocity space. As defined 
by Eqs. (3) and (4), Tq^^ and characterize temperature 
fluctuations along the line of sight and are similar to those 
introduced to explain the To([0 iii])/Tc(H i) discrepancy, 
which are given in Eqs. (1) and (2). In Eqs. (3) and (4), we 
can use 7^, (X*+), which is the intensity of a CEL used to 
determine Tc{v), as weights instead of /t,(II/3). However, if 
significant density variations are present along the line of 
sight, it may be problematic to use /«(X*+) since its density 
dependences differ in high- and low-density regions. 

Using Eqs. (3) and (4), we calculated Tq^^ and for 
three PNe, IC 2501, IC 4191, and NGC 2440. Our data 
were acquired using the Las Campanas Observatory (LCO) 
Baade 6.5m telescope with the MIKE echelle spectrograph, 
yielding a resolution of ~ 25,000 (see Sharpee et al. 120071 
for a detailed description of the observations). The results 
are given in Table [TJ For IC 4191 and NGC 2440 the values 
of Tq^v and t"^ are measured using the [O ill] AA4363, 5007 
lines. For the low-excitation PN IC 2501, the slit was placed 
on the edge of the compact O'^^ regions and the thermal 
broadening provides a dominant contribution to the profiles 
of the [O III] lines. For IC 2501 we obtained To,„ and tl 
using instead the single-ionization lines [N il] AA5754, 6678 
lines. For the calculations, the reddening coefficient and the 
electron densities presented by Sharpee et al. (2007) were 
adopted. As depicted in Fig. [3l the temperature variations 
in the velocity space are minor for the three PNe. 

Based on the values ofT^[0 iii]) and To(II i) pre- 
sented by Sharpee et al. p007p . we determined To and 
t^ using Peimbert's formulae (Peimbert il967j) . The calcu- 
lations were based on the assumption that the 0^+ and 
H+ regions are identical, which is appropriate for our stud- 
ied PNe since their 0+/(0+-|-0^+) abundance ratios are 
consistently lower than 0.1 (Sharpee et al. 2007). The re- 
sultant To and t^ are given in Table [TJ If the discrepancy 
between Te([0 iii]) and Te(H i) is caused by small tem- 
perature fluctuations, the t"^ measured by Te([0 ill]) and 
Te(H i) should be identical to those measured along the 
line of sight. However, Table [1] shows that the values of t^ 
are consistently lower than those of t^ and the values of To_„ 
are consistently higher than those of Tq, seemingly favoring 
the presence of extremely-cold nebular components where 
the [O III] lines cannot be excited but the H i recombination 
spectrum is significantly strengthened. 

This conclusion, however, needs to be treated with cau- 
tion. The value of t^, defined by Eq. (4) represents the tem- 
perature variations in velocity space. We should bear in 
mind that the validity of the argument that t^ is equal to 
the mean square temperature variation parameter along the 
line of sight requires two conditions to be satisfied: a) neb- 
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ular expansion-velocity dominates the broadening of line 
profiles; b) the expansion-velocity is monotonic along the 
nebular radius, as found in most of PNe. Otherwise, tf, pro- 
vides only a lower limit to the temperature fluctuations 
along the line of sight. In any case, the comparison of 
and t'^ provides an opportunity to investigate the presence 
of extremely-cold components. This paper is beginning such 
a project. Higher-resolution spectroscopic data of PNe with 
a high expansion velocity will be invaluable for this study. 

3.1.3. ORL profiles versus CEL profiles 

Figure also depicts that if temperature variations exist, 
the O II lines show a profile that differs from that of the 
[O III] lines, which can be attributed to the fact that CELs 
and ORLs weight high- and low-temperature regions, re- 
spectively. In the low velocity regions, the electron temper- 
ature is relatively high. As a result, the O ii lines have a 
deeper "gap" between the two peaks than the [O ill] lines. 
In the high velocity regions which correspond to the line 
wings, the electron temperature is slightly higher, depress- 
ing recombination lines. Consequently, the O ii lines are 
narrower than the [O ill] lines. The O ii line profiles differ 
most significantly from those of the [O ill] A4363 line than 
of the [O III] A5007 fine, since the [O ill] A5007 line has a 
lower excitation temperature than the [O ill] auroral line. 

High spectral-resolution observations of PN NGC 7009 
show that the [O ill] A4363 CEL is broader than the 
O II ORLs by a factor of about 1.5 (Barlow et al. 
I2006p . For their observations, a 0.9" x 0.9" image-slicer 
was placed on a bright-edge region located 5.6" north- 
west of the central star. We attempted to model the ob- 
served CEL/ORL width discrepancy. However, we find that 
chemically-homogeneous nebula models are unable to re- 
produce a CEL/ORL width ratio that is greater than 1.2 
because our modellings are unable to achieve such large 
temperature variations. The modelled CEL/ORL width ra- 
tio is essentially independent of the adopted velocity field 
and is insensitive to the assumed density distribution since 
the [O III] A4363 CEL and the O ii ORLs both have high 
critical densities. To match the observed CEL/ORL widths 
properly, the electron temperature is required to increase 
sharply outwards, which is difficult to achieve using pure 
photoionization modelling of gaseous nebulae. The most 
direct explanation for the observed large CEL/ORL width 
discrepancy is that the 0^+ CELs and ORLs originate in 
different nebular components which have separate kine- 
matic fields. Liu (200ir) proposed that evaporating plan- 
etesimals within the nebulae might produce cold H-deficient 
inclusions where ORLs are greatly enhanced and CELs are 
depressed. If this is the case, these H-deficient components 
should move in relatively-stable orbits, and thus have a 
lower velocity dispersion than the diffuse nebular gas. As 
an alternative, it remains possible that heating by shock 
waves in the outer regions may be significant and leads to 
a sharp temperature increase in the high- velocity regions. 
Furthermore, Stasinska & Szczerba (2001) showed that pho- 
toelectric heating by dust grains can cause large tempera- 
ture fluctuations in the presence of density inhomogeneities. 
However, in the case of a constant dust-to-gas ratio in the 
entire nebula, the temperature is more efficiently enhanced 
in the inner regions. If the CEL/ORL width discrepancy is 
caused by heating by dust, these dust grains must be pre- 
dominantly located in the outer zones of the ionized gas. 
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Fig. 4. Modelled results illustrating the line profiles for an 
outwardly-decreasing density distribution: a) the predicted 
line profiles of [S ii] A6716, [S ii] A4731, and [S ii] A4070; 
b) the assumed velocity field; c) the radial density distribu- 
tion of hydrogen atoms; d) the radial S+/S abundance-ratio 
distribution. 



3.2. Emission line profiles as a probe of density variations 
3.2.1. Modellings 

The determination of nebular density structure is critical 
to the investigation of ORL/CEL abundance discrepancies. 
Viegas & Clegg (|1994p pointed out that the presence of high 
density regions could lead to an overestimate of ^^^([O ill]) 
and consequently an underestimate of CEL abundances 
since the A5007 nebular line has a far lower critical den- 
sity than the [O ill] A4363 auroral line and is more signif- 
icantly suppressed by coUisional de-excitation in-high den- 
sity regions. However, Liu et al. (|2000p claimed that den- 
sity inhomogeneities can be ruled out as the cause of the 
ORL/CEL abundance discrepancies because no correlation 
between the abundance discrepancies and critical densities 
was found. Line profiles can be used to probe density vari- 
ations along the line of sight. 

The [S ii] AA6731,6716 doublet lines are classical 
electron-density diagnostics. We therefore simulate the pro- 
files of the [S ii] lines under the assumption that density 
decreases outwards and increases outwards, as shown in 
Figs. [4] and [H respectively. We assume that the ionizing 
star has a temperature of 30,000 K. For the nebulae that 
we model, is the dominant ionization-state of sulphur. 
Figures [4] and [5] indicate that the the profiles of the [S ii] 
AA6731,6716 doublet lines differ if density variations are 
present. The [S ii] A6731 line has a higher critical density of 
~ 4 X lO^cm-3 than the A6716 line (- 1 x 10^ cm^^), and 
hence is less-affected by collisional de-excitation in high- 
density regions. For a density that decreases and increases 
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Fig. 5. Same as|3]but for an outwardly- increasing density 
distribution. 



outwards, the [S ii] A6731 line is then stronger in low- and 
high- velocity regions, respectively. As shown in Figs. |4] and 
El the profile discrepancy is particularly remarkable when 
the [S ii] AA6731,6716 lines and the [S ii] A4076 hue are 
compared because the latter has a much higher critical den- 
sity (~ 10^ cm~'^). However, we note that temperature fluc- 
tuations may contribute partially to the profile discrepancy 
between these density diagnostic lines. 

Apart from for the density distribution, the basic as- 
sumptions are identical in each model. Figure [5] shows a 
larger separation between the two emission peaks than in 
Fig. 31 This is because in the case of outwardly-increasing 
density distribution, emission lines weight the high- velocity 
regions more highly. 

The analysis methods are applied to other density di- 
agnostics, such as [O II] AA3726,3729, [CI iii] AA5517,5537, 
and [O li] AA4711,4740. As an example. Fig. shows the 
profiles of the [Ar iv] AA4711, 4740 lines for high-excitation 
PNe with a stellar temperature of 150 000 K under the as- 
sumption of outwardly decreasing and increasing density 
distributions. In the two sorts of density structures, the 
relations between the two line profiles are completely the 
opposite. The [Ar iv] A4711 line has a lower critical den- 
sity than the [Ar iv] A4740 line, and thus weights lower- 
density regions more highly. Hence, if the nebula has a 
negative density gradient, the [Ar iv] A4711 line weights 
the outer regions more, where the expansion velocity is 
larger, which produces a broader profile in the [Ar iv] A4711 
line compared to the [Ar iv] A4740 line. For a similar rea- 
son, the [Ar iv] A4711 line has a narrower profile than the 
[Ar iv] A4740 line for an outwardly increasing density dis- 
tribution. A comparison of emission line profiles for various 
ionic species provides tighter constraints on nebular density 
structure. 



1.3 t-r 




Velocity (Itm/s) 

Fig. 6. The profiles of the [Ar iv] AA4711,4740 lines for 
outwardly-decreasing (upper panel) and increasing (lower 
panel) density distributions. 



3.2.2. Density variations in the velocity space 

Using the [S ii] AA6731, 6716 line profiles, we calculated the 
density distribution in velocity space, Nc{v), for the PNe 
IC 2501, IC 4191, and NGC 2440 (see Sharpee et al. lMJTI 
for the observation details). The results are shown in Fig. [71 
An inspection of the figure clearly indicates that the density 
variations in the velocity field are small (less than 1 dex). 
However, for the same reasons discussed in Sect. 13.1.21 for 
Tc{v), density variations in the velocity field may have been 
smoothed and only represent a lower limit to the nebular 
density variations unless line broadening is dominated by 
the expansion velocity and the velocity field is monotonic 
along the nebular radius. 

Figure |7| shows that the compact PN IC 2501 has an al- 
most homogeneous density structure with an average den- 
sity of ^ 10^ cm~^. The slight increase in density at 30 km/s 
is not real considering the large errors introduced by the 
weakness of the line wings. IC 4191 is a compact and irreg- 
ular PN. Figure [71 shows that it has an asymmetrical den- 
sity structure. A condensation is evident around 10 km/s 
on the red side, which has a density about twice as high as 
the other nebular regions. We cannot find a counterpart on 
the blue side. NGC 2440 clearly shows a negative density 
gradient across the nebula. The density at the outer edge 
of the nebula is lower by an order of magnitude than at the 
centre. 

We find that the [S il] A4076 line has a different profile 
from the [S ii] AA6731, 6716 lines for all three PNe. However, 
the [S ii] A4076 line was not taken into account in the cur- 
rent analysis because its upper level has a higher excitation 
energy than the [S ii] AA6731,6716, and hence we cannot 
separate the contribution of temperature variations to the 
difference between their profiles. 
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Fig. 7. Electron-density distribution in the velocity space 
of the PNe, IC 2501, IC 4191, and NGC 2440. 

3.3. Emission line profiles of IC 418 

To demonstrate how line profiles could be used to measure 
nebular physical conditions, we considered data for IC 418 
to complete a detailed study of line profiles by attempting 
to match our models to observational data. The analysis 
of a larger number of PNe will appear in a separate pa- 
per. High-resolution echelle spectra of the PN IC 418 were 
obtained in 2001 with the 4 m Blanco telescope at Cerro 
Tololo Inter- American Observatory (Sharpee et al. I2004p . 
covering a wavelength range from 3500 to 9865 A. IC 418 is 
a low-excitation PN [F([0 ill] A5007)/F(H/3) - 2] and has 
an approximatively-spherical shell. The observations were 
carried out using a ^ 1" slit width, yielding a spectral reso- 
lution of ~ 30, 000. The spectrograph slit was aligned north- 
south. The central star was placed along a line perpendic- 
ular to the spectrograph slit so that the slit centre was 
roughly midway between the central star and outer edges 
of the nebula. A large number of ORL and CEL profiles are 
well resolved. 

Hyung et al. p994p constructed a photoionization 
model of IC 418 that can interpret the UV data fairly well. 
For our model, we initially took the basic parameters used 
by Hyung et al. 119941 (distance, properties of the central 
star, chemical composition, nebular geometrical structure, 
etc.), which were then adjusted slightly to achieve an op- 
timal fit to optical observations. The density distribution 
of IC 418 is a topic of ongoing debate. Hyung et al. (|1994p 



presented that the nebula may consist of two shells, an in- 
ner high-density one and an outer lower-density shell. A 
conflicting conclusion was obtained by Gesicki et al. (1996]) 
who found that the density increases smoothly outwards by 
fitting the surface-brightness profiles for H/3, [O ill] A5007, 
and [N ii] A6584. Our model assumed a single shell with 
a homogeneous density of A'h = lO^cm"^. Below we show 
that this issue can be settled by comparing the profiles of 
density diagnostic lines. 

Once the photoionization model provides satisfactory 
fits to the observed integrated line intensities, we calcu- 
late the line profiles by assuming a velocity field and com- 
pare the predicted results with the observed ones. For the 
comparison, we used a number of strong lines from differ- 
ent ionic species, including H i. He i, [S ii], [N ii], [O i], 
[O ii] , [O III] , and [Ne ill] . The observed and predicted line 
profiles are presented by Fig. [8] Low-ionization lines have 
generally broader profiles than high-ionization lines, which 
implies that the expansion velocity is increasing outwards. 
To reproduce the observed line profiles, the acceleration is 
required to increase in the outer regions, as shown in Fig. [51 
For the fitting, no turbulence broadening is needed. Figure[5] 
shows that our model can account for most of these strong 
lines. 

Figure [9] plots the FWHMs of the lines as a function of 
ionization energy, as predicted by our model and observed 
by Sharpee et al. (|2004p (see their Fig. 7). Good agreement 
between the predictions and the observations is achieved. 
The velocity field produced is in general agreement with 
the results of Gesicki et al. (|1996p . but has a larger outer 
velocity. Given that Gesicki et al. (|1996p did not consider 
the [O l] line, their measured outer velocity is unlikely to be 
reliable. The sharply-increasing velocity in the outer regions 
is associated with the shock at the ionization front, which 
is consistent with the predictions of hydrodynamic models 
(e.g. Perinotto et al. 11^55)) . 

The [O ii] A7330 line can be explained by our model, 
while the observed [O ii] A3726 line is broader than the pre- 
dicted one. This is unlikely to be due to temperature vari- 
ations since the [N ii] AA5754, 6584 doublet lines is repro- 
duced well. The discrepancy, therefore, can be attributed to 
density variations within the PN. The [O il] A3726 line has 
a critical density about three times lower than that of the 
[O ii] A7330 line, and thus originates in low-density regions. 
It follows that the outer regions have a lower density than 
the inner regions, in contrast to the results of Gesicki et al. 

mm- 

For the [O i] A6300 line, our model reproduces the pro- 
file in the high-velocity range well. But in the low-velocity 
range, the observed flux is much higher than predicted. 
Morisset & Stasinska l|2006ap also encountered difficulty in 
fitting the low-velocity range of the [O i] A6300 line pro- 
file. A possible explanation is that neutral gas is plenti- 
fully abundant within the ionized regions. This conjecture 
was supported by Williams et al. (|2003p . who compared 
absorption- and emission-line spectra of IC 418 and found 
evidence that absorbing gas may be present. If we enhance 
the [O i] emissivity in the inner regions by a factor of ~1.5, 
an excellent match between the predicted profile and ob- 
served one is achieved, as shown the Fig. [S] 

For observations of low-excitation PN, the slit was 
placed in the outer edge of the O ii regions where the pro- 
jected velocity along the line of sight is quite low. Therefore, 
the narrow high-ionization lines, such as [O ill] and [Ne ill]. 
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Fig. 9. FWHM as a function of the ionic species' ionization 
potential for the PN IC 418. The solid line with triangles 
is the prediction of our model. The box denotes the zone 
occupied by the observations of Sharpee et al. (,2004,) . 




Fig. 8. Observed (solid lines) and predicted (dotted lines) 
strong emission lines of IC 418. The assumed velocity field 
is shown in the lower panel. For the [O i] A6300 line, the 
dashed line represents the predicted profile by refined emis- 
sion distribution, as discussed in the text. 



are dominated by thermal and instrumental broadening. 
The [O III] AA4363, 5007 lines are therefore unsuitable for 
deriving temperature variations in velocity space. 

We compare the predicted ORL profiles with the obser- 
vations. Figure [To] displays a few examples for such a com- 
parison. Generally, our model shows good agreement with 
observations of the high-photoionization lines, C ii, N ii, 
O II, and N ill. IC 418 has a small ORL/CEL abundance 
discrepancy [(02+/H+)orl/(02+/H+)cel - 1.3; Sharpee 
et al. 12004] , and is therefore an unsuitable target to inves- 
tigate the materials in which ORLs originate. 

Figure [10] demonstrates that our model provides a rea- 
sonable match to the O i A5513, but is unable to explain the 
profiles of O i A7773 and N i A8223, whose widths appear 
to be broader that those of the predicted lines. As pointed 
out by Sharpee et al. (|2004p . fluorescence excitation may 
contribute substantially to the intensity of these lines, and 
the radiation of the lines may originate predominantly in 
the outer neutral regions. 



Fig. 10. Same as Fig. [S] but for recombination lines. 
4. Conclusions 

This paper addresses the potential of emission line pro- 
files to probe the physical conditions of PNe. Our results 
show that if temperature or density variations are present 
within the nebulae, lines with different excitation tempera- 
tures or critical densities might show significantly different 
profiles even though they originate in similar ionic species. 
We show that line profiles provide a new way to investi- 
gate the CEL/ORL abundance-discrepancy problem. For 
this purpose, we require high-resolution and deep spectro- 
scopic data. In the ideal case, the line width contributed 
by instrumental broadening should be smaller than that by 
the thermal broadening of heavy-element lines. This techni- 
cal requirement corresponds to a spectral resolution higher 
than 75, 000. 

We present the electron temperatures and densities in 
the velocity space of a sample of PNe. No significant tem- 
perature or density variations are found. This provides a 
lower limit to the temperature and density inhomogeneities 
along the line of sight. 

We attempt to reproduce the observational data of two 
PNe, N GC 6153 and NGC 7009, acquired by Barlow et al. 
pOOe]) . which show that the [O ill] A5007, the [O ill] A4363, 
and the O ii recombination lines have different profiles. 
We find that a pure photoionization model of chemically- 
homogeneous gas can explain the [O ill] A5007/A4363 pro- 
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file discrepancies, but cannot explain tlie [O iii]/0 ii pro- 
file discrepancies. We thus conclude that CEL and ORL 
may originate in kinematically-difFerent nebular compo- 
nents. Alternatively, the [O iii]/0 ii profile discrepancies 
could be caused by extra heating in the outer regions. In 
this work, ID computations were completed. Line profiles 
may significantly depend on nebular geometrical structures, 
as proposed by Morisset & Stasinska (2006a). However, for 
lines from the same ionic species, nebular geometrical struc- 
tures provide an effect on the profiles of a similar degree, 
and thus hardly affect our results. 

We construct line profiles of the approximatively- 
spherical PN, IC 418. Our model can explain most of the 
observed line profiles. A velocity field that sharply increases 
outwards is revealed. We find that the [S li] and [O ii] den- 
sity diagnostic lines have very different profiles, which imply 
that large density variations are present within the PN. Our 
model shows that in the low-velocity regions of the [O i] 
A6300 line, the predicted flux is lower than the observed 
one, which indicates the existence of neutral clumps within 
the ionized regions. Generally, the profiles of CELs and 
ORLs are in good agreement. However, IC 418 is a young 
PN which has a small CEL/ORL abundance-discrepancy, 
and thus is unsuitable for studying the abundance prob- 
lem. For further studies, high-quality spectroscopic observa- 
tions of PNe with large CEL/ORL abundance-discrepancies 
should be invaluable. 
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